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Abstract 

A quantum fractal is a wavefunction with a real and an imaginary 
part continuous everywhere, but differ entiable nowhere. This lack of 
differentiability has been used as an argument to deny the general va- 
lidity of Bohmian mechanics (and other trajectory-based approaches) 
in providing a complete interpretation of quantum mechanics. Here, 
this assertion is overcome by means of a formal extension of Bohmian 
mechanics based on a limiting approach. Within this novel formu- 
lation, the particle dynamics is always satisfactorily described by a 
well defined equation of motion. In particular, in the case of guidance 
under quantum fractals, the corresponding trajectories will also be 
fractal. 



1 Introduction 

Quantum mechanics is the most powerful theory developed up to now to 
describe the physical world. However, its standard formulation, based on 
statistical grounds, does not provide an intuitive insight of microscopic phe- 
nomena as classical mechanics does for macroscopic ones. For example, the 
evolution of a system cannot be followed in the configuration space by means 
of well denned, individual trajectories. On the contrary, the wavefunction 
associated to such a system extends to the whole available space, describing 
the probability for the system to be located at each (space) point at a certain 
time. 

In order to obtain a more intuitive picture of quantum phenomena, al- 
ternative approaches relying on the concept of trajectory have been proposed 
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PP . One of them is Bohmian mechanics El El • This theory is not merely a 
reinterpretation of the standard quantum mechanics despite its equivalence 
at a predictive level, but a generalization of classical mechanics that accounts 
for quantum phenomena. Hence, since Bohmian mechanics formally rests on 
the same conceptual grounds as classical mechanics, the notion of trajectory 
(causally describing the particle evolution) can also be applied in a natu- 
ral way to study of microscopic world without contradicting the statistical 
postulates of the standard quantum mechanics. 

Recently, it has been argued |3] that trajectory-based theories, like Bohmian 
mechanics, fail in providing a complete interpretation of quantum mechanics. 
In particular, these theories could not satisfactorily deal with wavefunctions 
displaying fractal features, the so-called quantum fractals [HI 13 E] • The ex- 
istence of this class of wavefunctions in quantum mechanics has important 
consequences from a fundamental viewpoint: despite the coarse-graining re- 
strictions implied by Heisenberg's uncertainty principle, these wavefunctions 
constitute the proof that fractal objects can appear in quantum mechanics 
as well as in classical mechanics. Indeed, Wocik et al. jS] pointed out that 
quantum fractals could be experimentally constructed by considering heavy 
atoms or ions in macroscopic traps, where the number of energy levels would 
be large enough to observe scaling properties at least up to several orders of 
magnitude (what is considered a physical fractal). On the other hand, Ama- 
natidis et al. have observed theoretically this fractal behaviour during 
the ballistic and diffusive evolution of wave packets moving in tight-binding 
lattices, a study of interest in quantum information theory and quantum 
computation. 

Here it is shown that the incompatibility between Bohmian mechanics 
and the existence of quantum fractals can be easily avoided by reformulating 
the particle equation of motion. This reformulation is based on a limiting 
approach, where the wavefunction is expanded in a series of eigenvectors of 
the Hamiltonian. If the wavefunction is not fractal, but regular, the tra- 
jectories are the same as directly computed by means of standard Bohmian 
mechanics. However, if the wavefunction presents fractal features, its non- 
differentiability forbids a direct calculation of the trajectories, which can be 
obtained, on the other hand, by means of the limiting approach proposed 
here. Thus, this formal extension of Bohmian mechanics, based on a novel 
reformulation of the particle equation of motion, provides a causal picture 
for any arbitrary wavefunction, regular or fractal. 

The organization of this paper is as follows. In order to make the paper 
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self-contained, a survey on quantum fractals is given in section |21 The funda- 
mentals of Bohmian mechanics and its generalization to deal with quantum 
fractals are presented in sectional The application of the new concepts intro- 
duced in this work is illustrated in detail in section0]by means of the problem 
of a non-relativistic, spin-less particle of mass m in a one-dimensional box. 
This simple, integrable problem can be considered a paradigm of fractals 
appearing under conditions not necessarily related to a chaotic dynamics [8] . 
In section El the question of the unbounded energy for quantum fractals and 
its interpretation in terms of quantum trajectories is discussed. Finally, the 
main conclusions derived from this work are summarized in section |3 

2 Quantum fractals 

A general method [Hj to construct quantum fractals with an arbitrary fractal 
dimension consists in using the quantum analog of the Weierstrass function 

oo 

W(x) = ^b r sm(a r x), a > 1 > b > 0, ab > 1 (1) 

r=0 

the paradigm of continuous fractal function. Thus, for example, in the prob- 
lem of a particle in a one-dimensional box of length L (with < x < L), 
solutions of the Schrodinger equation can be constructed as 

R 

$ t (x; R) = A R n r(s ~ 2) smip^x/H)^^ (2) 

r=0 

with 2 > s > and n > 2. Here, p n>r = n r nh/L and E n , r = p^ r /2m 
are, respectively, the quantized momentum and the eigenvalue associated 
to the eigenvector that corresponds to the quantum number n' = n r ; and 
Ar is the normalization constant. The wavefunction (J2J) is continuous and 
differentiate everywhere, however, the one resulting from the limit 

$t(x) = lim $ t (x;R) (3) 

is a fractal object JT] in both space and time. 

This method to generate quantum fractals basically consists (given s) 
in choosing a quantum number, say n, and then considering the series that 
contains its powers, n' = n r . There is another alternative (and related) 
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method [U| to obtain quantum fractals based on the presence of discontinuities 
in the wavefunction. In this case, although the initial wavefunction can be 
relatively regular, fractal features emerge due to the perturbation that the 
discontinuities cause on the wavefunction along its propagation. 

An illustrative example of this kind of generating process is a wavefunc- 
tion initially uniform along a certain interval, I = x 2 — X\ < L, inside the 
box mentioned above, 

*o(s) = | ^' . (4 ) 

[ 0, elsewhere 

The Fourier decomposition of this wavefunction is 

2 °° 1 

V (x) = — - - [cos(p n xi/H) - cos(p n x 2 /h)] sm(p n x/fr) (5) 

and its time-evolved form is 
2 00 1 

%(x) = — -= - [cos(p n x 1 /h) - cos(p n x 2 /h)] sm(p n x/h) e~ iEnt/h . (6) 
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n=l 



As can be noticed, this wavefunction is equivalent to assume r = R = 1 
in (J2J), and sum over n, from 1 to N, obtaining the quantum fractal in 
the limit iV — > 00. This equivalence is based on the fact that the Fourier 
decomposition of \l/o gives precisely its expansion in terms of the eigenvectors 
of the Hamiltonian in the problem of a particle in a box. However, this is not 
general, since the Fourier decomposition and the expansion of \l/o m a basis 
of eigenvectors of the Hamiltonian are not equivalent when V is not constant 
along x. 

The fractality of wavefunctions like (jHJ) or (jUJ) can be analytically esti- 
mated [H] by taking advantage of a result from Fourier analysis. Given an 
arbitrary function 

K 

f(x) = J2^~ iKX (7) 

re=l 

its real and imaginary part are fractals (and also \f{x)\ 2 ) with dimension 
Df — (5 — P)/2 if its power spectrum asymptotically (i.e., for K — > 00) 
behaves as 

kr-N^ ( 8 ) 
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with 1 < (3 < 3. Alternatively, the fractality of f(x) can also be calculated 
by measuring the length, C, of its real or imaginary part (or \f{x)\ 2 ) as a 
function of the number of terms, K, considered in the generating series (J7J). 
Asymptotically, the relation between L and K is given by 

C(K) oc K ?- 1 (9) 

which diverges for f(x) being a fractal object. Notice that increasing the 
number of terms that contribute to f(x) is analogous to measuring its length 
with more precision, since its structure is gradually better determined. 

A remarkable feature that characterizes quantum fractals is that the ex- 
pected value of the energy, (H), of these wavefunctions is unbounded. This 
is related to the fact that the familiar expression of the Schrodinger equation 

ihdt%(x) = H%(x) (10) 

does not hold in general as happens when *&t( x ) is a quantum fractal. In 
this case, neither the l.h.s. of equation (|1U|) nor its r.h.s. belong to the Hilbert 
space. Hence, the equality is not formally correct, and the applicability of 
this equation fails. On the contrary, since each term of the series satisfies 
this equation, the identity 

[H-ihd t ]%(x) = (11) 

which also represents the Schrodinger equation, still remains valid. When 
this happens, ^(x) is called jS] a solution of the Schrodinger equation in a 
"weak" sense. 

3 Quantum fractal trajectories 

The fundamental equations of Bohmian mechanics are commonly derived by 
writing the system wavefunction in polar form 

%( x ) = p] /2 {x)e lSt{x)/h (12) 

with p t = \^t\ 2 being the probability density and St the (real-valued) phase, 
and substituting it into the Schrodinger equation (fTUj). This leads to two 
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(real-valued) couple differential equations 



l+V.f^U. (13) 



dt \ m 

f + M£ + v + ft -a 

Equation (|T3|) is a continuity equation that ensures the conservation of the 
flux of quantum particles. On the other hand, equation (j!4j). more inter- 
esting from a dynamical viewpoint, is a quantum Hamilton- Jacobi equation 
governing the motion of particles under the action of a total effective poten- 
tial V t eS = V + Qt- The last term in the l.h.s. of this equation is the so-called 
quantum potential 

h 2 vV /2 

*=-£;^p- (15) 

This context-dependent, non-local potential determines together with V the 
total force acting on the system. 

In the classical Hamilton-Jacobi theory, St represents the action of the 
system at a time t, and the trajectories describing the evolution of the system 
correspond to the paths perpendicular to the constant-action surfaces at 
each time. Analogously, since the Schrodinger equation can be rewritten 
in terms of the Hamilton-Jacobi equation (|14j) . S t can be interpreted as a 
quantum action satisfying similar mathematical requirements as its classical 
homologous. The classical concept of trajectory emerges then in Bohmian 
mechanics in a natural way, defining the particle trajectory as 

it = Y^i = 1 Im r*-iv* t l . (i6) 

mm 

Since in Bohmian mechanics the system consists of a wave and a particle, it 
is not necessary to specify the initial momentum for the particles, as hap- 
pens in classical mechanics, but only their initial position, xq, and the initial 
configuration of the wavefunction, ^q. The initial momentum field is pre- 
determined by \l/o via equation (JHJ), and the statistical predictions of the 
standard quantum mechanics are reproduced by considering an ensemble of 
(non-interacting [12] ) particles distributed according to the initial probabil- 
ity density, p . 

Equation (fTB^l is well defined provided that the wavefunction is continu- 
ous and differentiable. However, this is not the case for quantum fractals. 
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This is the reason why one might infer a priori that Bohmian mechanics 
is an incomplete theory of quantum motion [5] unable to offer a trajectory 
picture for this type of wavef unctions. This apparent incompleteness can 
be nevertheless "bridged" by taking into account the decomposition of the 
quantum fractal as a sum of (differentiable) eigenvectors of the corresponding 
Hamiltonian, and then redefining equation in a convenient way. 

Since regular wavefunctions are particular cases of quantum fractals for 
which the fractal and topological dimensions coincide, the new, generalized 
equation of motion will be applicable to any arbitrary wavefunction, \l/ t . Such 
a wavefunction can be expressed as 



N 



y t (x;N)=J2cnUx>- iE " t/h (17) 

n=l 

with N — > oo, and where £ n (x) is an eigenvector of the Hamiltonian with 
eigenvalue E n ; in the case where the wavefunction is constituted by a limited 
number M of eigenvectors, c n = for n > M. Accordingly, the quantum 
trajectories evolving under the guidance of (fT7|) are defined as 

x t = lim Xjsr(t) (18) 



with xjv(t) being the solution of the equation of motion 



XN{t) = — lm 
m 



(19) 



Observe that this reformulation of Bohmian mechanics is not totally 
equivalent to the conventional one. The calculation of trajectories is not 
based on S t , which cannot be trivially decomposed, in general, in a series of 
analytic, differentiable functions, as happens with \l/ t . Thus, the existence 
of trajectories is directly postulated taking into account equations (|T5jl and 
(119)) rather than equation (|14|). For regular wavefunctions both formulations 
coincide due to the differentiability of St. Whereas, when dealing with quan- 
tum fractals, the particle equation of motion is only well defined within this 
reformulation, and gives rise to quantum fractal (QF) trajectories. The frac- 
tal dimension of these trajectories can be determined by means of equation 
(JHJ), now C referring to the QF-trajectory length. 
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4 A numerical example 



The problem of a highly delocalized particle inside a one-dimensional box 
(i.e., a particle with the same probability to be found everywhere inside the 
box) illustrates fairly well the concepts described in the previous sections. 
The wavefunction representing the state of such a particle is given by equation 
with X\ — and X2 = L. Taking this into account, equation © becomes 

%(x) = _JL e - iE i*A - sm(p n x/h) e"^ 1 * (20) 
n vL nodd n 

where u ny i = (E n —Ei)/h. In the numerical calculations, L = m = h = 1 (in 
arbitrary units, a.u.). 

The probability density, p t , associated to the wavefunction (|271j) is a pe- 
riodic function of time, with period T = 2tt/uj^^ = mL 2 /2iTh. To show that 
this is the periodicity of p t is relatively easy. At t = T, the arguments of the 
interference terms contained in p t are io n ^T = 2n(n 2 — l)/8, with n > 1. In- 
deed, since n is always an odd integer, it can be written as n(k) = 2(k — 1)+3, 
with k > 1, and then u) n (k),iT — k(k + l)ir. This result shows that at t — T 
any argument is always an integer multiple of 2ir, and therefore the minimum 
time elapsed between two consecutive recurrences is precisely T = 2ti/uj 3 i 

(* = !)• 

Despite the periodicity of p t , the wavefunction (j2U|) is not truly periodic 
due to the common time-dependent phase, <p t = —Eit/H, multiplying the 
sum (for example, the wavefunction undergoes a delay of —tt/4 after each 
period). This delaying phase is a general feature for any wavefunction ex- 
pressible as dHJ), but has no consequences from a quantum trajectory view- 
point. Equation (|18p is invariant under space-independent factors added to 
the phase St, since 

x N [S' t ] = x N [S t ] (21) 

when S' t = St + s(t). Here, in particular, s(t) = (ft- The invariance of the 
quantum motion with respect to such factors is consistent with the fact that 
two wavefunctions that differ in a phase factor represent the same state in 
standard quantum mechanics. From now on S t will refer to the phase of (|20j) 
without the factor s(t). 

The profiles along x of p t and S t are displayed, respectively, in figures 
da) and (b) at two different times. These functions display a fractal shape 
or a revival (characterized by a step-ladder shape) depending on whether the 
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Figure 1: Probability density (a) and phase (b) associated to a highly delo- 
calized particle in a box at t — (thin solid line) and t = 0.7 T (thick 

solid line), (c) Measure of the fractal dimension of the probability densi- 
ties displayed in part (a). To compare, measures of the fractal dimension 
of initial probability densities associated to triangular (T) and parabolic (P) 
wavefunctions are also shown. 
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time is an irrational or a rational fraction of the period, respectively. The 
fractal-revival alternation manifests the Cantor-set structure [TO] of (|2(Jj) 
along time (i.e., its real and imaginary part display an infinite number of 
alternating fractal and revival profiles along time). 

The revivals are characterized by the well-known Gibbs phenomenon re- 
lated to the Fourier decomposition of discontinuous functions, which does 
not affect the quantum motion. Apart from this, as seen in figure [Ha), these 
revivals also present regions close to the boundaries of the box where pt van- 
ishes at certain times; the most dramatic case happens at t = T/2, when 
Pt{x) 7^ only in the interval 0.5 < x < 0.75. These nodal regions are very 
important from a dynamical viewpoint. Since St is not well defined in these 
regions (observe that S t is not represented for x < 0.05 and x > 0.95 in figure 
[TJb)), particles avoid them. 

The fractal nature of \l/ t is quantified by applying equation to p t . 
The logarithm of the length (C) of pt as a function of the logarithm of N is 
represented in figure ^c) for the two cases considered in part (a). As clearly 
seen, log 10 £ is proportional to log 10 N in the fractal case, resulting a fractal 
dimension Df = 1.49, which is in excellent agreement with that obtained by 
Berry [HJ using equation (jHJ). On the other hand, as expected, the length 
corresponding to the revival approaches a constant saturation value. The 
eventual growth observed in the graph is related to the slow convergence of 
p t to a step-ladder structure. If other regular probability densities with no 
discontinuities are considered, the convergence is much faster. This happens, 
for example, when one considers that is a triangle (T) or a parabola (P), 
both centred at x c = 0.5. In these cases, also represented in figure life), 
the saturation is reached relatively faster, since only few eigenvectors are 
necessary to obtain an excellent convergence. The slower convergence in the 

(T) (T) 

case of Pq is due to the non-differentiability of \&q at x c , which implies a 
larger number of eigenvectors in the sum. 

The complex space-time structure generated by p t along its evolution, 
the so-called fractal quantum carpet [IS], can be easily understood by means 
of the QF-trajectories, which provide a causal description for such a pattern. 
As seen in figure I21^a), these trajectories manifest the symmetries displayed 
by ^t, the guiding wave. Thus, in the case of the reflection symmetry with 
respect to x c , the trajectories started at one side of the box (to the left 
or right of x c ) do not ever cross to the other side. This effect due to the 
single-valuedness of St, which avoids the trajectories to cross at the same 
time, can be compared with a hard-wall scattering problem; an ensemble 
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Figure 2: (a) QF-trajectories associated to a highly delocalized particle in 
a box. (b) Measure of the fractal dimension of a sample of QF-trajectories 
with initial positions: xq = 0.01 (■), xq = 0.1 (•), xq = 0.4 (A), xq = 0.49 
(□), x = 0.499 (o), and x = 0.5 (A). 



of particles initially moving towards the wall will display similar features to 
those observed in figure EJa) (see, for example, JEHU)]). Here, the particles 
cannot cross the point x c , acting like a fictitious wall, and therefore they 
bounce backwards describing trajectories symmetric with respect to t = T/2. 
This inversion of the particle momentum is related to the second kind of 
symmetry that affects the wavefunction propagation: the change of sign of 
St during the second half of the period. Moreover, unlike classical trajectories 
(and also caused by the single-valuedness of S t ), not all QF-trajectories can 
reach the wall, but will move parallel to it. This is a nice manifestation of 
the effects caused by the quantum pressure under fractal conditions. 

In figure El^b), the logarithm of the length of several QF-trajectories, 
log 10 £, is given as a function of log 10 N. As clearly seen, the converge to 
proportionality is faster for those QF-trajectories started at intermediate 
positions, between the boundary and the centre of the box. Notice also 
that, since the trajectory started at x c is not a fractal, its length does not 
depend on N. Independently of the initial position (and with the exception 
of the trajectory started at Xq = 0.5), the fractal dimension of any trajectory 
asymptotically approaches the same value, -D/ — 1.50, which coincides with 
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that found for p t . 

5 Causal considerations about the infinite- 

ness of (H) 

As seen in section |21 the expected value of the energy becomes infinite 
for quantum fractals. In the wavefunction (J2J), for example, the condition 
2 > s > gives rise to a divergent series in (H) when the limit Q is 
taken into account. This property, unavoidable when dealing with quantum 
fractals, is related to their infinite scaling behaviour [S], which is lacking 
in (everywhere and anytime) regular wavefunctions. In the case of revivals 
of wavefunctions with discontinuities, like ()20|) . which are regular only at 
certain times, (H) also remains unbounded because an infinite number of 
eigenvectors is necessary to recreate the discontinuities. In this way, the dis- 
continuities can be understood jB] as perturbations that propagate along the 
box in time, leading to the formation of the quantum fractals. 

A more physical insight on the infiniteness of (H) can be gained by ap- 
pealing to the trajectory formulation introduced in section El In Bohmian 
mechanics, the particle energy is given by 



Except in the case of particles associated to eigenvectors of the Hamiltonian, 
E t does not conserve in time in general, although the average energy of an 
ensemble of particles, initially distributed according to po, 



does it, in agreement with the standard quantum mechanics. The classi- 
cal analog of two couple particles can help to easily understand this fact; 
although the energy of each particle varies along time due to a continuous 
transfer between both, the total energy will remain constant. Therefore, 
though Bohmian particles are independent > the presence of the quantum 
potential in equation (JHj) leads to sort of non-local coupling or dependence 
between each particle and the rest from the ensemble (whose evolution is 
described by equation (H2J)). 




(22) 




(23) 
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In the example used in section 0] the potential energy is V = at any 
time. Thus, the initial total energy of almost all particles is zero, since they 
are distributed according to a constant probability density (Qq = 0) and the 
wavefunction is real (So = 0). Only at the boundaries of the box Q = oo 
due to the discontinuity. This infinite amount of energy is stored up in the 
particles located at x = e and x = L — e with e — ► + , which constitute 
energy reservoirs. These energy reservoirs are the set of particles located at 
the discontinuities of whenever a revival emerges, and not only at t = 0. 

The dynamical evolution of the system described by the wavefunction 
(120}) can be explained in terms of this initial non-homogeneous energy dis- 
tribution among particles as a function of their initial position. Although 
the particle distribution is homogeneous in space, the particles are not in 
quantum equilibrium, but subjected to an infinite gradient of energy at the 
boundaries of the box. This gradient leads to a strong, symmetric energy 
flow going from the boundaries towards the centre of the box that makes 
particles to move as shown in figure |2f a). 

The relationship between the time dependence of the energy and the QF- 
trajectory dynamics can be better understood, without loss of generality, by 
recalling a simpler example consisting in assuming N = 3 in the wavefunc- 
tion (|20|). An ensemble of trajectories illustrating the dynamics associated 
with this case is shown in figure Efa); these trajectories can be regarded as 
coarse-graining envelopes of the QF-trajectories shown above. Although the 
trajectories have not been distributed according to p , they provide an insight 
on how p t evolves, which is represented in figure HHe) at three different times: 
t\ = 0.15 (dotted line), t 2 = 0.35 (dashed line), and £3 = 0.47 (thick solid 
line), in units of T. The total and kinetic energy, and the quantum potential 
are represented, respectively, in figures HDT)-(h) at these three times. On the 
other hand, these magnitudes are also respectively displayed as a function of 
time in figures Efb)-(d) for three different trajectories with initial conditions: 
Xq = 0.01 (thick solid line), x Q 2 ' = 0.05 (dashed line), and xf^ = 0.3 (dot- 
ted line). As seen, the total energy can be, alternatively, positive or negative 
along time since it does not conserve. 

The quantum potential acting on a particle depends on the structure 
of p t at each position of the particle. Thus, local minima in p t translate 
into negative wells in Q t that particles avoid [TH], moving towards regions 
with positive values of Q t , or, at least, presenting local maxima. The most 
dramatic case occurs when p t has a node (see pt at £3), manifested as a 
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Figure 3: (a) Quantum trajectories associated to the wavefunction PU)) with 
N = 3. The total and kinetic energy, and the quantum potential are repre- 
sented, respectively, in parts (b), (c) and (d) for three different trajectories 
of panel (a) with initial conditions: x^ = 0.01 (thick solid line), x$ = 0.05 
(dashed line) and x = 0.3 (dotted line), (e) Probability density at three 
different times: t\ = 0.15 (dotted line), t 2 = 0.35 (dashed line) and = 0.47 
(thick solid line), in units of T. These times are indicated in panel (a) by 
parallel, horizontal lines. The total and kinetic energy, and the quantum 
potential at these three times are represented, respectively, in parts (f), (g) 
and (h). 
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singularity in Q t . The intense forces around these regions make particles to 
move extremely fast apart from them ^7], provoking peaks in their kinetic 
energies, as seen in figure H2g). This behaviour is not observed, however, 
for those wells that appear in the central region. As commented above, this 
is because the particle cannot cross the point x c in these cases, but keeps 
moving for a certain time close to it until the quantum pressure decreases 
sufficiently, and can move backwards. 

In this way, at t± very few particles remain close to the borders of the box 
(the maxima of Qt are relatively narrow), and most of them move towards 
the maxima located around x ~ 0.3 and x ~ 0.7, as can be seen in figureEfa). 
This gives rise to the two important peaks in p t ; see figure Efe). At t 2 , the 
marginal maxima occupy a wider extension, however they are relatively high 
in energy, and therefore most of the particle flow directs towards the central 
maxima. As a consequence, pt displays an important central maximum, and 
two secondary, marginal maxima; see figure 0(e). Finally, at £3, as seen in 
figureEfe), most of the particles are collected in the centre of the box, between 
x ~ 0.25 and x ~ 0.75, impelled by the strong forces around x ~ 0.16 and 
x ~ 0.84, respectively. Moreover, since Q t reaches high values at the borders, 
very few particles will remain in such regions. 

According to this analysis, it is clear that the three kinetic energy curves 
shown in figure 0(g) display peaks on the minima of Q t , as would happen in 
a classical situation. However, the transient trapping observed along x c has 
a purely quantum nature, since there is no physical (classical) potential that 
may contribute to it [IH[T5]. By following the sequence tx-t^-t^ one can see 
that Q t progressively increases at the borders and develops deep wells that 
confine the particles within the central part of the box. In other words, the 
quantum pressure increases from the borders of the box towards the centre, 
pushing the particles towards x c , and obliging them to move along this point 
for some time (approximately, half a period). 

Taking into account the ideas exposed above, the analysis of a single 
particle dynamics becomes much simpler. For example, as seen in figure 
0(a), the trajectory x[ is initially slightly pushed away, towards x c , by Q t 
until it reaches a turning point, and then moves backwards. Because of this, 
two peaks are observed in its kinetic energy, the second smaller than the first 
because the particle does not turn back to the original position. The turning 
point, like in classical mechanics, is characterized by a zero value of K t . The 
fact that Q t reaches its minimum at the turning point can be understood 
as an appearance of a non-crossing wall (similar to that at x c ) avoiding 
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the particle to go beyond it. On the other hand, between t/T = 0.3 and 
t/T = 0.7, the particle remains on top of the plateau seen in figure IH^h), and 
it is almost at rest (there is only a very slight oscillation at about t/T = 0.5). 

(2) 

The same is analysis applicable to the trajectory x t , although the changes 
in its velocity are much more relevant, mainly at about t/T = 0.5, when the 
particle undergoes the strong force due to the singularity in Q t . Finally, for 
the trajectory x\ the double peak is only observed at half of its evolution 
unlike the two previous cases. This is because this particle does not reach any 
turning point during the first part of its evolution, but only a sudden force 
that pushes it towards x c in a fast manner. Once in the trapping region (with 
the highest quantum pressure) the particle oscillates, and finally undergoes 
another sharp forcing that separates it from the neighbourhood of x c . 

In the light of the previous analysis, one can conclude that the variation in 
time of the energy can be understood as a regulating mechanism that adjusts 
the particle motion in such a way that it results consistent with the evolution 
of the wavefunction. The discussion applied to the trajectories guided by a 
three-state wavefunction also remains valid in the case of QF-trajectories. 
However, the motion adjustment takes place in a relatively faster manner, 
since particles will reach an infinite amount of turning points along their time 
evolution. Therefore, Q t and K t will display very deep wells and very sharp 
peaks, respectively, and the total (average) energy, given by equation (J2HJ), 
will diverge. 

6 Conclusions 

The consistent picture of quantum motion provided by Bohmian mechanics 
relies on a translation of the physics contained within the Schrodinger equa- 
tion into a classical-like theory of motion. This transformation from one 
theory to the other is based on the regularity or differentiability of wave- 
functions. Therefore, it does not hold for quantum fractals, non-regular 
solutions of the Schrodinger equation. A priori, this seems to be a fail- 
ure of Bohmian mechanics in providing a complete explanation of quantum 
phenomena, since quantum fractals would not have a trajectory-based rep- 
resentation within its framework However, taking into account the fact 
that Bohmian mechanics is formally equivalent to the standard quantum 
mechanics, this incompleteness results quite "suspicious". 

By carefully studying the nature of quantum fractals, one can under- 
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stand the source of such an incompatibility. These wavefunctions obey the 
Schrodinger equation in a weak sense jS], i.e., given the wavefunction as a 
linear superposition of eigenvectors of the Hamiltonian, the Schrodinger equa- 
tion is satisfied by each eigenvector, but not by the wavefunction as a whole. 
This is because the eigenvectors are always continuous and different iable 
everywhere, unlike quantum fractals, which are continuous everywhere, but 
differentiable nowhere. Taking this into account, a convenient way to express 
any arbitrary wavefunction, regular or fractal, is in terms of a superposition 
of eigenvectors of the Hamiltonian. This procedure is particularly important 
in those circumstances where the differentiability of the wavefunction is going 
to be invoked, like in the formulation of trajectory-based quantum theories 
like Bohmian mechanics. 

In order to have a truly consistent particle equation of motion, Bohmian 
mechanics must be then reformulated in terms of an eigenvector decompo- 
sition of the wavefunction instead of considering the latter as a whole (as 
happens in standard Bohmian mechanics). The resulting generalized equa- 
tion of motion, defined by a (convergent) limiting process, is valid for any 
arbitrary wavefunction, and provides the correct Bohmian trajectories. In 
the case of quantum fractals, one obtains the desired trajectory-based pic- 
ture at the corresponding limit. Whereas, if the wavefunction is regular, the 
trajectories determined by means of this procedure will coincide with those 
given by the standard Bohmian equation of motion. This novel generaliza- 
tion thus proves the formal and physical completeness of Bohmian mechanics 
as a trajectory-based approach to quantum mechanics. 

The trajectories associated to quantum fractals are also fractal. This 
explains both the formation of fractal quantum carpets and the unbounded 
expected value of the energy for quantum fractals. Although the example of a 
particle in a box has been used here to illustrate the peculiarities of quantum 
fractals, the analysis can be straightforwardly extended to continuum states 
[2] or other trajectory-based approaches to quantum mechanics, like Nelson's 
theory of quantum Brownian motion Moreover, this kind of analysis can 
be of practical interest in the study of properties related to realistic systems, 
like those suggested by Wocik et al. [S] and Amanatidis et al. 0, providing 
moreover a causal insight on their physics. 
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